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In theories with large extra dimensions, supernova (SN) cores are powerful sources of Kaluza-Klein 
(KK) gravitons. A large fraction of these massive particles are gravitationally retained by the newly 
born neutron star (NS). The subsequent slow KK decays produce potentially observable 7 rays and 
heat the NS. We here show that the back-absorption of the gravitationally trapped KK gravitons 
does not significantly change our previous limits. We calculate the graviton emission rate in a 
nuclear medium by combining the low-energy classical bremsstrahlung rate with detailed-balancing 
arguments. This approach reproduces the previous thermal emission rate, but it is much simpler 
and allows for a calculation of the absorption rate by a trivial phase-space transformation. We 
derive systematically the dependence of the SN and NS limits on the number of extra dimensions. 
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I. INTRODUCTION 

Theories with large extra dimensions are a recent 
alternative to solve the hierarchy problem of particle 
physics Within a certain class of models, 

the most restrictive limits on the size of the extra dimen- 
sions derive from the supernova (SN) emission of Kaluza- 
Klein (KK) gravitons, particles with an essentially con- 
tinuous spectrum of masses that are a generic feature of 
the new theory. Even though these new particles interact 
very weakly, i.e. the strength of ordinary gravitons, the 
number of modes and thus the size of the extra dimen- 
sions is constrained by the requirement that SN 1987A 
did not emit more KK gravitons than is compatible with 
the observed neutrino signal duration (2,7 A 9i lOj- 

In the simplest models, KK gravitons are stable ex- 
cept for their slow, gravitational-strength decay into pho- 
tons, neutrinos, and other standard particles. There- 
fore, the decays of KK gravitons produced in all cos- 
mic SNe will contribute to the measured cosmic 7-ray 
background, providing more restrictive limits than the 
SN 1987A energy-loss argument 

Long after the parent SN has exploded, a neutron 
star (NS) will continue to shine in 7-rays because a 
large fraction of the KK gravitons remains gravitation- 
ally trapped — most of them were produced with masses 
near the kinematical production threshold and thus with 
small velocities. Therefore, a NS is embedded in a halo 
of KK gravitons that shines in 100 MeV 7-rays. The 
EGRET non-observation of such signatures from nearby 
NSs thus provides restrictive limits The strongest 

constraints yet, however, are from avoiding excess heat- 
ing of certain old NSs which otherwise could not cool to 
their observed low surface temperatures T^]. 

These arguments depend on the assumption that KK 
gravitons are not re-absorbed by inverse nuclear brems- 
strahlung, a question that we failed to address when we 



raised our argument |l2l| . A crude estimate of the re- 
absorption effect gives a time scale larger than the ages 
of the stars that we used for our limits, but close enough 
to warrant a more detailed calculation. Moreover, if the 
re-absorption time scale is shorter than the KK decay 
lifetime, NSs will be heated directly by KK absorption 
rather than absorbing the decay products. The excess 
heat would increase and our bounds would improve. 

The first goal of the present note is to calculate KK 
graviton absorption in a NS. Essentially this is done by 
phase-space transforming the emission rate. This exer- 
cise is particularly simple and transparent if one writes 
the graviton emission rate in a form that separates the 
response of the thermal nuclear medium from the phase 
space of the radiation. This approach significantly sim- 
plifies the original calculation of the emission rate of 
Hanhart, Phillips, Reddy and Savage and illuminates 
the nature of the approximations made. To achieve the 
same level of precision we only need Weinberg's classi- 
cal bremsstrahlung rate |l3j | together with the principle 
of detailed balancing. This approach closely follows pre- 
vious calculations of neutrino or axion emission from a 
SN core based on the general properties of the thermal 
medium's response functions [l4L [l5l|. Therefore, while 
the main goal of our derivation is to obtain the graviton 
absorption rate, its derivation is an illuminating exercise 
in its own right. 

Our second goal is to extend our limits to the general 
case of n extra dimensions. Previously, the limits on the 
compactification scale were explicitly stated only for the 
n < 3 cases. We use this opportunity to show how the 
various SN and NS limits scale with n. 

In Sec. II we derive the emission and absorption of 
ordinary gravitons in the low-energy limit from a nuclear 
medium. In Sec. Ill we extend these results to the case 
of KK gravitons. In Sec. IV we revisit the SN and NS 
constraints, and in Sec. V we summarize our results. 
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II. NUCLEAR GRAVI-BREMSSTRAHLUNG 
A. Differential Energy-Loss Rate 

The dominant graviton emission process from a SN 
core is nuclcon-nucleon bremsstrahlung NN NNg 
where for now g stands for an ordinary graviton. In the 
limit of soft radiation the bremsstrahlung rate is deter- 
mined by the measured nuclear scattering cross section 
alone, i.e. the nuclear scattering event and the associated 
bremsstrahlung process factorize . In the soft limit the 
details of the scattering process do not matter, only the 
sudden change of the nucleons' energy- momentum tensor 
is responsible for the radiation. 

Assuming that two non-relativistic nucleons collide 
with the initial CM momenta Pi = Pi = — P2 to the 
final states p/ = P3 = — P4, the differential amount of 
emitted gravi-bremsstrahlung energy is [l3| 



duj 



16 Gn pi 



5n 



2(P.-P/)' 



M2 



(1) 



where Gn is Newton's constant, M the nucleon mass, and 
uj the graviton energy. We may also write this in terms of 
the total initial and final kinetic nucleon energies, Ei = 
pI/2M + pI/2M = pf/M and Ef = pl/2M + pl/2M = 
p'j/M, and in terms of the CM scattering angle as 



dE„ 



16 G^ 



E 



Ej - 2E,Ef cos^ GcM 



Energy conservation implies Ef = Ei — lu so that 



dEa 



) Gpv 



{Ei + Ef)^sm^e 



CM 



(2) 



(3) 



div 5n 

an approximation used in Ref. 0. 

If the emission of soft gravitons {uj 0) is viewed as 
a classical process, i.e. the colliding nucleons as external 
sources, energy conservation does not apply, Ef = Ei, 
and we get the usual flat bremsstrahlung spectrum. 



dE^ 
duj 



32 G^ 



Ef sin^ e 



CM 



(4) 



Augmenting this classical result with detailed balancing 
will give us the correct emission rate up to 0{lu) as in 
the more complicated treatment of Ref. • 

Next, we introduce the nucleon scattering amplitude 
A which is normalized such that the differential cross 
section for neutron-neutron scattering is 



da_ _ M^TT^ \A\^ 
dfl " (27r)6 



(5) 



Note that \A\'^ includes a spin sum over initial- and final- 
state neutrons while da/dil is averaged over initial and 
summed over final spins. If one takes \A\'^ to be indepen- 
dent of scattering angles, the total cross section is 



apja^ 

327r 



25 mb. 



(6) 



A factor ^ for identical final-state nucleons was included 
and the numerical value was quoted in Ref. as a good 
approximation for the conditions of interest. 

With these ingredients the differential energy-loss rate 
of a neutron medium in the soft limit is 



dQ 
doj 



= j dTS\A[' 



32 G^ 



Ef sin^ e 



CM 



(7) 



where S* = i is a statistics factor for identical particles in 
the initial and final state while dV symbolizes the neutron 
phase space integration, including all thermal occupation 
numbers, blocking factors and the energy-momentum 5 
function. True to our soft-radiation approximation, the 
graviton energy-momentum does not appear in this 5 
function. Without the term in square brackets Eq. {Tj) is 
the rate of nn collisions per unit volume. 

We now read off the dependence of the energy-loss rate 
Q on various input parameters. If the medium is non- 
degenerate, the nn collision rate per unit volume is the 
neutron number density squared times the cross section 
times an average thermal relative velocity, the latter be- 
ing proportional to (T/Af)^/^ with T the temperature. 
Further, from the bremsstrahlung rate we have an aver- 
age Ef which is proportional to T^. A typical emitted 
graviton energy is of order T so that finally 



cx 



GNcrn?,r^/2M-V2 



(8) 



where is the density of baryons (here neutrons). 

We will perform explicit calculations only for the case 
of non-relativistic non-degenerate neutrons. In that case 
the phase-space integral can be transformed to one over 
CM momenta 0. Integrating a quantity F over the nu- 
cleon phase space yields 



dTF^n^[-\ 



X / dui I duf(uiUf) I e~^' S{ui ~ Uf) 
Jo Jo 

-1 
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X 



dcos8cM F 



(9) 



where Ui j — Ei j/T. We have explicitly kept the energy 
d function for later convenience. 

Applying this result first to the nn scattering rate per 
unit volume we insert F — S\A\'^. Assuming that the 
scattering amplitude is independent of the CM energy 
and the scattering angle we find explicitly 



1/2 



"■B 
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Wrol 



(10) 



where the factor i accounts for identical initial-state par- 
ticles. The corresponding final-state factor is already in- 
cluded in the total cross section. 
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Next we perform the analogous integral for the differ- 
ential energy-loss rate and find 



dQ 
duj 



r 



128 2 

nn — — L^N-t 
OTT 



256 Gnct^rTS/^ 



(11) 



57r3/2 

for the "classical" bremsstrahlung losses of the medium. 

B. Total Rate from Detailed Balancing 

The differential rate derived thus far is not yet use- 
ful for calculating the total Q. To this end we observe 
that the emission, absorption, or scattering of any radi- 
ation by any medium is described by a suitable dynam- 
ical structure function. This is a quantum correlator of 
those medium operators which couple to the relevant ra- 
diation. For example, for neutrino pair emission, pair 
absorption, or scattering, we need the nucleon spin cor- 
relator because neutrinos couple primarily to the nucleon 
spin 0, 01 . For graviton emission we need a dynami- 
cal structure function describing the microscopic fluctu- 
ations of the medium's energy-momentum tensor. 

However, in the soft limit when the energy transfer uj 
is small we may use the long-wavelength limit where the 
radiation's momentum transfer is ignored. The energy- 
loss rate of any type of radiation is then of the form 



2w(27r) 



■ UJ S{—uj) 



(12) 



where S{uj) is a function of the energy transfer alone. 
Apart from overall coefficients, S is the dynamical struc- 
ture function in the long- wavelength limit. The integral 
is over the invariant phase-space of the radiation, one 
factor CO accounts for the energy carried by the radiation. 
Note that from the medium's perspective a negative lu 
is energy lost, a positive uj energy gained. Therefore, 
bremsstrahlung emission involves S{—uj). 

The only property of relevance to our present discus- 
sion is the detailed-balancing requirement 



S{- 



S{uj) e 



ojT 



(13) 



which is a general consequence of the non-commuting na- 
ture of the correlated quantum operators describing the 
thermal medium, i.e. it strictly is a quantum effect. This 
condition ensures that the efficiency of emission and ab- 
sorption of radiation is such that in equilibrium the ra- 
diation will reach a thermal distribution function. As an 
immediate consequence the most general structure func- 
tion is of the form 



S{uj) = 5(w) ■ 



(14) 



where Siio) = \ \S(lS) -|- S'(— w)] is symmetric in uj. 



From Eq. lO we may now write the differential 
energy-loss rate in the form 



dQ_ 

duj 



Sico) 



(15) 



(27r)2 l + e'^/T' 
Comparing with Eq. (|ll|l we recognize that we may write 



(16) 



where s(x) is a dimensionless even function normalized 
to s(0) — 1 and 



Sn — 



1024 ttV^ G^an\T^I^ 
5 Ml/2 

1.547 X 10^^ erg cm"^ s"^ MeV"^ ^30^3 • 



(17) 

Here, T30 = T/BO MeV and pa = p/3 x Vi>^^ g cm"^. 
Since S(lS) is even apart from the "trivial" detailed- 
balancing factor, we obtain the full differential energy- 
loss rate up to first order in a; = LojT even though we 
have calculated the bremsstrahlung rate only to zeroth 
order in LojT . 

Our approach of using the concept of the structure 
function kills two birds with one stone. We obtain the 
differential emission rate up to order w/T on the basis of 
the strictly classical soft-radiation bremsstrahlung rate. 
Second, the inverse process, graviton absorption, is sim- 
ply a different phase-space integral over the same struc- 
ture function and thus trivial to extract. 

A calculation up to 0{lo/T) is the best one can do in 
the framework of the soft-radiation approximation. A 
more accurate calculation would require modeling the 
nucleon interaction in detail because the scattering and 
bremsstrahlung parts of the process no longer factorize. 
Of course, a more precise calculation would also require 
taking many-body effects into account, i.e. one would 
need to calculate properly the dynamical structure func- 
tion, an impossible task at present. 

In the absence of any more precise insights concerning 
the behavior of the structure function we may use the 
simplest approximation compatible with our level of ap- 
proximation and take s{x) = 1, dropping all higher-order 
terms. In this case we find for the total emission rate 



Q 



512 In 2 

57r3/2 



(18) 



The coefficient is numerically 12.75. 

The corresponding results of Ref. are recovered if 
in the nuclear phase-space integral Eq. Q we include 
the graviton energy, i.e. energy conservation now reads 
5(ui — Uf — x) with X = Lo jT . Moreover, we must include 
the expression Eq. (|3J| for dEg/duj. One first performs 
J dhj to remove the 5 function and then integrates the 
remaining expression. The coefficient of Q is then found 
to be 47104/(525 7r3/2) sa 16.11, somewhat larger than 
our result. 
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We can also extract the structure function implied by 
this treatment. To this end we leave the J duj integra- 
tion open and rather do J dui to remove the 5 function. 
Collecting all factors we find 



s{x) 



1 + e- 
48 



1+48^ 



duf{u 



(19) 



However, s{x) is not analytic, the coefficient of the x'^ 
term diverges. This expression indeed agrees with ours 
up to 0{x). Put another way, up to 0{uj) Ref. ex- 
plicitly recovers the detailed-balancing condition that we 
used as input information. 

We can go through the same exercise for dEg/duj of 
Eq. © and find s{x) = 1 + (5/16) a;^ + 0(a;^). This 
again deviates from the other results only at 0{x^). 

In summary, even though one has used the soft- 
radiation approximation which is valid only for uj <^ T, 
one obtains an estimate of the emission rate which likely 
is correct within a numerical factor of order 1. Still, this 
approach has essentially the status of an "educated di- 
mensional analysis." A real calculation would require 
knowledge of the relevant dynamical structure function 
of the nuclear medium, not just its long- wavelength, soft- 
radiation limit. 



C. Inverse Bremsstrahlung 

We are now in a position to calculate the graviton ab- 
sorption rate by a nuclear medium. For simplicity we 
may start with the energy- loss rate Eq. (|12|l . leave out 
the graviton phase space integration except for the fac- 
tor l/2a;, remove the factor lo which represents the emit- 
ted energy, include a factor 1/2 because the energy-loss 
rate is summed over two graviton polarizations, and sub- 
stitute to — + ~Lu in the structure function because the 
medium now gains energy. Then we find 



rg(^) 



4a; 

_ 256 ttV^ Gn<jp^ r'/' : 

5 l + e 

= 3.7 X 10"^^ s-^pI Tg/^c^r^ 



(20) 



where = p/3 x lO^^ g cm^^, T30 = T/SO MeV, and 
u!iao — w/100 MeV. We have assumed s{x) = 1 and in 
the numerical expression also e~^^^ = 0. 



III. KALUZA-KLEIN GRAVITONS 

A. Large Extra Dimensions 

Next we assume that space has n large extra dimen- 
sions in which gravity can propagate. String theory sug- 
gests a total of 11 space-time dimensions so that there 



are at most 7 large extra dimensions in addition to our 
4-dimensional space-time. Therefore, we will limit our 
explicit calculations and results io \ <n <1 . 

We further assume that the extra dimensions are 
toroidally compactified with an equal radius K for all 
of them. Put another way, the linear size of each ex- 
tra dimension is lixK and the volume of the compactified 
space is Vn = (27ri?)". 

Kaluza-Klein (KK) gravitons that propagate in the ex- 
tra dimensions with a momentum p will appear to us as 
having a mass m — p. The modes are discrete with a 
density of states Ki/(27r)" = R". A summation over all 
modes in the extra dimensions corresponds to summing 
over all masses of KK gravitons in our world. For large 
extra dimensions the modes are narrowly spaced so that 
the sum over masses transforms to an integral 



E 

modes 



dm 1 



where 



27r"/2 

r(«/2) 



(21) 



(22) 



is the surface of the n-dimensional unit sphere. 

Following the conventions of Ref. 6] we note that 
the fundamental scale of the new theory is the effective 
Planck scale M4+„ of the higher dimensional space. It 



relates to the conventional Planck scale Mp\ 
1.22 X IQi^ GeV by 



G 



-1/2 



M 



PI 



KM 



■2+n 
4+n 



(23) 



where A/pi — Mpi/^/Sn. Numerically we find 



M 



4+n 



1 TeV 



(5.922 X 10 



3Q\ — 



3.140 X 10"^" m 
R 



_ (24) 
for the relationship between M4+„ and R. 

Note that M4+„ is different from the energy scale 
M4+„ that was used by Hanhart et al. l9l.lir| and that we 
called M = M4+„ in our previous papers|lll IT^ . This 
energy scale was defined from the relation 



47rG'N = 



1 



RM 



1 

M2 



(25) 



With Eq. ^ this implies 

M = 2i/("+2) (27r)"/("+2) Mi+n . (26) 

This relationship allows one to translate our new limits 
on M4+„ into limits on the previous parameter M. In 
the second row of Table IVII we give explicit values for 
Af/Af4_|_„, i.e. one should multiply our new limits with 
this factor to obtain the corresponding limits on the old 
parameter M. 
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B. Absorption and Emission of KK Gravitons 



where 



For KK gravitons, the absorption rate will be similar 
to that for ordinary gravitons of the same energy. There- 
fore, we may use Eq. (pn|l directly for this purpose. 

For the energy-loss rate we must take account of the 
large number of available modes. For a fixed KK mode 
with mass m one has to include in Eq. H12I) a factor |^ 



19 11 m2 2 
18 ^ y ^ ^ 9 "w^ 



(27) 



that already includes a summation over polarization 
states. In addition, J d^k — > An Juikduj with k = 
{lJ^ — w^Yl"^ . The total emission rate then obtains by 
summing over all masses. With Eq. (|21|l one thus finds 



^0 



where we have defined 



, X VI. 2n1/2 /19 11 2 2 A 



(28) 



(29) 



using /i ~ m/oj. Performing the integral over masses 
explicitly yields 



9nR" 

(27r)2 J, 



dww"+2S'(-w) 



(30) 



where 



gn = ^n / duGn-lifJ.) 



5n— 19 -I- 18n + 



48 



r(^) 



(31) 

Note that go = 19/18 does not represent the emission of 
ordinary gravitons because of differences in the summa- 
tion over polarization states. 

With these results it is straightforward to state the 
differential energy-loss rate. We read from Eq. H3l)|l 



dQn So 



doj (27r)2 
Likewise 

dQn _ So 
dm (27r) 



s{uj/T) . (32) 



^n^iRTT (f) " F(m/T), (33) 



where 



Fiy) 



dx 



^-s{x)Go(^-) (34) 



with y = m/T and x = u>/T. For s{x) = 1 a simple 
approximation is F{y) « 2y°-^^ e~^'^'^^y that is good to 
about ±7% for O.K /i < 20. 

If we finally use our expression for the structure func- 
tion with the simplest assumption s{x) = 1 we find 



(35) 



= 2.35 X 10^^ erg cm^^ pi . (36) 

Here, T30 = T/30 MeV and ps = p/3 x IQi"^ g cm'^ and 

/„ = (l-2-")r(n + l)C(n + l), (37) 

valid for n> I. 

In Table U we give numerical values for [gnfn)'^^^ be- 
cause this is the effective n-dependent coefficient of R in 
the emission-rate formula. A simple approximation is 



(gnfn) 



1/r, 



(38) 



good to better than 1% in the range 1 < n < 10. 

The corresponding emission rate of Ref. |^, their 
Eq. (56), is systematically larger by factors of up to a 
few because their s{x) grows quadratically with x, in- 
creasing the emission of higher-energy gravitons. Recall 
that for the emission rate one needs the structure func- 
tion for Lu/T of a few while the soft-radiation approxi- 
mation is only valid for lu <^ T. Therefore, our result 
agrees with Ref. Q up to the claimed accuracy of either 
result. The soft-radiation calculation can not predict the 
emission rate to better than a factor of a few. 

It is straightforward to calculate the average values of 
several parameters of the emitted gravitons. In Table UTI 
we give explicit numerical values for the average energy 
(w) , the average mass (m) , and the average velocity (v) . 
We also give the retention fraction /trap in a typical neu- 
tron star and the average mass of the trapped gravitons 
as explained in the following Section. 

TABLE I: Various coefficients described in the text. 



n 
1 
2 
3 
4 
5 
6 
7 



(<7n/n)^/" 

1.7943 
2.4544 
2.8981 
3.2597 
3.5771 
3.8662 
4.1348 



0.3763 
0.9731 
1.1607 
1.2001 
1.1891 
1.1605 
1.1263 



TABLE II: Average properties of the emitted KK gravitons. 



n 


{lo}/T 


(m)/T 


(v) 


/trap 


(m)trap/r 


1 


1.1866 


0.5784 


0.8061 


0.0963 


1.0902 


2 


2.1923 


1.4054 


0.7052 


0.1816 


2.0176 


3 


3.1514 


2.2663 


0.6384 


0.2607 


2.9052 


4 


4.1060 


3.1516 


0.5892 


0.3342 


3.7917 


5 


5.0691 


4.0601 


0.5506 


0.4020 


4.6890 


6 


6.0429 


4.9888 


0.5191 


0.4643 


5.5993 


7 


7.0257 


5.9335 


0.4927 


0.5211 


6.5208 
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C. Neutron-Star Retention of KK Gravitons 

The KK gravitons are produced with barely relativistic 
velocities. Therefore, the SN core will retain a large frac- 
tion of them gravitationally after thermal production by 
the nuclear medium. In order to estimate the retention 
fraction we take the NS to be a homogeneous sphere with 
a density 3 x 10^^ g cm~^ and a mass A/ns — 1-4 Mq, cor- 
responding to a radius of Rns — 13.06 km and a surface 
gravitational potential of 



-0.1139. 



(39) 



We here use Newtonian physics rather than general rela- 
tivity. Our simple NS model leads to a radial dependence 
of the gravitational potential within the NS of 



U{r) = (Jns 



(40) 



where r is a dimensionless radial coordinate that varies 
from to 1 in the NS. A KK graviton produced at ra- 
dius r with energy lu remains trapped if its kinetic plus 
gravitational energy is negative, U{r) + < or 

w > m > mtrap(r) =uj[l + U{r)] , (41) 

where we have used that \U\ <C 1. 

The graviton number emissivity is given by Eq. H30|l if 
we drop one power of uj under the integral. The number 
emissivity of those KK gravitons that remain trapped, 
averaged over the entire NS, is found by replacing gn 
with g^^^ in these expressions. The latter is found if in 
Eq. we substitute 



I dfi ...-^3 [ drr^ [ 

Jo Jo J u. 



(42) 



where /^trap = Wtrap/^^ = 1 + U{r). 

For fixed n the distribution of emitted masses is a 
universal function of the parameter fi = m/uj, indepen- 
dently of the energy oj. The velocity of a KK graviton 
is V ~ {1 — l^'^y^^ and thus only a function of fi so that 
the distribution of KK velocities is also independent of 
their energy. Therefore, the trapping fraction is indepen- 
dent of the KK energy, in every energy class u the same 
fraction is retained given by gn^^/gn- Explicit trapping 
fractions for different values of n are shown in Table ITll 

As another consequence the average energy of the 
trapped KK gravitons is identical with that of the es- 
caping ones. However, the average mass of those that 
remain trapped is much larger than the overall average. 
We give explicit results for (m)trap in Table ^ These 
values are only slightly smaller than (w), reflecting the 
fact that the trapped particles are slow-moving and thus 
have only small amounts of kinetic energy. 



IV. CONSTRAINTS ON LARGE EXTRA 
DIMENSIONS 

A. Supernova 1987A 

We begin with the classical SN 1987A energy-loss 
limit on the emission of KK gravitons 0, S Il0| . 
The duration of the SN 1987A neutrino signal pre- 
cludes that too much energy was emitted in an in- 
visible channel. For several cases involving axions or 
right-handed neutrinos, self-consistent cooling calcula- 
tions were performed to determine the relationship be- 
tween the allowed coupling strength of the exotic parti- 
cles and the duration of the neutrino burst. The result 
can be summarized by the simple criterion that the ex- 
otic energy-loss rate of a nuclear medium at a density 
of 3 X lO^"' g cm^'^ and a temperature of 30 MeV should 
not exceed 10^^ erg g~^ s~^ [ig.[l7| . This criterion is not 
crude or arbitrary, but is calibrated by numerical simu- 
lations of different authors with different codes and for 
different cases and reproduces these results surprisingly 
well. The numerical studies of Ref. for the specific 
case of KK gravitons confirms once more the accuracy of 
the simple criterion. 

We have only considered a non-degenerate neutron 
medium. For the realistic case of a mixture of pro- 
tons and neutrons the emission rate would be somewhat 
larger, while degeneracy effects would slightly decrease it. 
The soft-radiation calculation is only accurate to within 
a factor of a few. Therefore, it is not worth worrying 
about the exact infiuence of the mild nucleon degener- 
acy or the chemical composition. The simple energy-loss 
criterion together with Eq. (|35() thus implies 



R < 6.58 X 10~^^ m (1.28 x 10^ 



17\1/" 



(gnfn) 



-1/n 



(43) 



We give the explicit limits for 1 < n < 7 in Tabled The 
corresponding limits on Af 4+„ are given in Table EH 

For n = 2 our limit R ;$ 0.96 fj,m compares with 
0.70 iim of Ref. and 0.66 /xm of Ref. For n = 3 
our limit i? ^ 1.14 nm compares with 0.83 nm of Ref. 
and 0.8 nm of Ref. !l^. The small differences are per- 
fectly in line with the different approximation made in 
the calculation of the emission rate and with our using 
the simple energy-loss criterion rather than a specific nu- 
merical simulation as in Ref. jlOj . 



B. All Cosmic Supernovae 

In the toroidal compactification scheme assumed in our 
discussion the only decay channels for KK gravitons are 
into those standard-model particles that are kinemati- 
cally allowed, notably into 27, e~e~^ and vD. The cou- 
pling is of gravitational strength so that the mean lifetime 
is very large 

1 ^ „„ /100MeV\^ ,,,, 
^27 = T^^e-e^ = ~ 6 X 10^ yr . 44 

2 \ m I 
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This is comparable to the age of the universe for masses 
in the 100 MeV range that would be typical for those KK 
gravitons emitted by a SN core. 

The photons from these decays from all past SNe would 
contribute to the cosmic 7-ray background as measured 
by the EGRET instrument 19], constraining the amount 
of KK gravitons that may have been emitted by all cosmic 
SNe 11]. The measured diffuse 7 background falls with 
energy roughly as E~^. We have explained in Ref. 
that this particular power-law behavior implies that the 
limit on the KK emission by SNe is nearly independent of 
the assumed core temperature and nearly independent of 
the KK mass spectrum as long as the masses are so large 
that most of the gravitons would have decayed within a 
Hubble time. This is the case for T > 20 MeV and n>2. 

In Ref. 11] we found that the EGRET data imply that 
a typical SN core must not lose more than about 0.5% of 
its energy in KK gravitons. With a realistic SN rate the 
limit would be more restrictive by factors between 10 and 
100. In the spirit of deriving conservative limits it is thus 
justified to assume that for all n not more than 1% of 
the energy loss is allowed that we used for the SN 1987A 
limit. Put another way, the SN 1987A limits on R scale 
with (10~^)^/". The corresponding explicit limits for 1 < 
n < 7 are given in Tables Ivl and IVII respectively. 



photons per graviton. We also assume that for every 
decaying KK graviton we have = u>/2, i.e. we ignore 
the narrow energy distribution of the decay photons from 
non-relativistic gravitons. Collecting all factors we pre- 
dict an EGRET flux at Earth of 



(27r)2 



47rd2 



X 2 / dujS{-uj) 

'2Eo 



(imw"+'^G„+3(m/a)) 
(47) 



where Eq = 100 MeV. The / dm integral is explicitly 
w"+4/i„ with 

5n(278 + 223n + 48n2-)-3n3)0Fr(f) , ^ 
hr, = — ^ , — — . (48) 



384r(2±li) 



Therefore, we have 



where 



i{u) ^ dx 



1 



(49) 



(50) 



C. Supernova Remnant Cas A 

The SN remnant Cas A probably corresponds to Flam- 
steed's SN of 1680, but in any case is so young that the 
cloud of emitted KK gravitons would still appear as a 
point source to EGRET, even ignoring the fact that a 
large fraction of them is gravitationally retained [T]| . 
The absence of an EGRET source at the location of 
Cas A implies that the flux of decay photons at Earth 
from this source is limited by 



and 



£;^>100 McV 



< 10"^ cm"2 g-i 



(45) 



In order to predict the photon flux at Earth we begin 
with the KK emissivity Eq. (|28(l . To get the number flux 
we must include a factor under the integral. Further, 
we must include the decay rate r2-y, involving the time- 
dilation factor m/u} and t^^ from Eq. (|^ . i.e. 



$n = 



1 FNsAiNs745oT3 



(27r)2 47rd2 
= 7.1 X 10-24 cm-^s-iTg^'ps. 



(51) 



We given explicit values for (f7„/i„)^/" in Table 

Comparing with Eq. (|45|l then gives us the limits on 
R and M4+„ shown in Tables Ivl and IVII respectively, for 
T = 30 MeV. In contrast to the previous results they now 
depend rather sensitively on the SN core temperature. 
Reducing T to 20 MeV degrades the _R-limit for n = 2 by 
about a factor of 10, for n = 7 still by about a factor of 2. 
The Mi^n limits are somewhat less sensitive. Changing 
from 30 to 20 MeV degrades the n = 2 limit by about 
a factor of 3.5, the n — 7 limit by about a factor of 2. 
Overall the Cas A limits are comparable to those from 
all past SNe. 



F 



27 = 7 



(46) 



D. Gamma Radiation from Neutron Stars 



with 7 = 5.3 X 10-24 s-i MeV-3. We further need 
to multiply with the volume Vns of the NS and with 
the time-scale At^s of KK emission, i.e. the cooling 
time scale of the nascent neutron star. We will use 
^NS = M-^s/p, assume a NS mass of 1.4 Mq, and take 
AiNS — 7.5 s. This time scale corresponds to the 
SN 1987A energy-loss argument, assuming that at most 
half of the total energy goes into KK gravitons. Finally 
we need to divide by iird"^ with d = 3.4 kpc the dis- 
tance to Cas A. We further need a factor 2 for two decay 



A large fraction of the KK gravitons emitted by a SN 
core are gravitationally retained so that every NS would 
be embedded in a halo of these particles. Therefore, NSs 
would be bright sources of 100 MeV 7-rays visible to 
EGRET. The flux of decay photons expected from a NS 
can be written in a way which is very similar to Eq. H47|l . 
except that the integral over dm is changed so that only 
the trapped gravitons are counted 



(52) 
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where 



8.1 X 10 



-23 



^30 Psrf; 



2 

kpc ■ 



(53) 



The function is defined as 



In{u) 



1 

r^dr 



l-;7(r) 



„ra+2 



1 + 



(54) 







In Table Iml we give /„ for n = 1-7 and T = 30 MeV. 

Comparing this flux prediction with the EGRET point- 
source limit Eq. H45|l for various old NSs then allows one 
to derive limits on R. In our previous paper 12] we 
used several old NSs for this argument. The most re- 
strictive limits obtained from the nearest neutron star 
RX J185634-3754 that was taken to be at a distance of 
60 pc [23 and from PSR J0953+0755 at 120 pc 
However, the distance of RX J185634-3754 has recently 
been adjusted upwards to 117 ± 12 pc [23|, so that the 
distance is now exactly the same as for PSR J0953-I-0755. 
This in turn means that the derived limits on R from the 
gamma-ray luminosity is exactly the same for the two 
stars. The limits on R and M4+„ thus derived are tabu- 
lated in Tables IVl and IVll respectively. 

The flux prediction of Eq. scales with the SN core 
temperature as T^^/^, just as for the Cas A argument. 
Therefore, the limits are equally sensitive to T as dis- 
cussed in the previous section. 



E. Neutron-Star Excess Heat 

1. Decay of KK gravitons 

Even more stringent constraints obtain from consider- 
ing the heating of old NSs b y th e trapped cloud of KK 
gravitons surrounding them |l2l |. When a graviton de- 
cays outside the NS, the fraction of the decay photons 
hitting the surface of the star is i-^/d^ — R'^g/d, where 
d is the distance of the graviton from the center of the 
star, and i?NS is the NS radius. 

We proceed in the following way to calculate the rms 
distance of a particle in the graviton cloud from the cen- 
ter of the star (d^)^/^. First we calculate the average 
total energy E = ^v^ + U{r) for the trapped gravitons. 
We then assume for simplicity that all gravitons are on 
purely radial orbits. Given the level of approximations 
entering this calculation this is certainly justified. By 
solving the equation of motion for a graviton with en- 
ergy i? on a radial orbit we get [d"^)^/"^ , and from this we 
can calculate the fraction, of the total KK-produced 
photon flux absorbed by the NS. In doing this we assume 
that all gravitons decay when they are at a distance of 
exactly (d^)^/^. Numerical values for (d^)^/^ and T are 
shown in Table Hvl for T = 30 MeV. 



TABLE III: Values of the parameter J„ 
(j>n{2Eo/T) and In{2EQ/T), taken for T = 
100 MeV. 



and the functions 
30 MeV and Eq = 



n 






I„{2Eo/T) 


J n 


1 


0.605 




7.56 X 10"^ 


0.292 


2 


4.933 




0.0443 


1.061 


3 


41.42 




0.274 


4.731 


4 


360.2 




1.797 


25.05 


5 


3.27 X 10^ 




12.48 


153.8 


6 


3.11 X 10** 




92.17 


1074.5 


7 


3.13 X 10^ 




727.3 


8428.2 


TABLE IV: Values of {d^)^/^ 


, tinsidc, and the reabsorption 


rate r,i. 










n 




T 


^inside 




1 


1.405 


0.149 


0.269 


1.22 X lO"^*^ 


2 


1.372 


0.158 


0.279 


2.06 X lO-^'^ 


3 


1.352 


0.164 


0.287 


7.18 X 10"^* 


4 


1.320 


0.174 


0.297 


3.44 X 10"^* 


5 


1.291 


0.184 


0.308 


1.94 X 10"^* 


6 


1.262 


0.195 


0.319 


1.21 X 10"^* 


7 


1.236 


0.206 


0.330 


8.14 X lO^^'^ 



Having calculated the fraction of the total flux ab- 
sorbed by the star, we can calculate the total energy per 
unit time absorbed by the star, L, as 



L = ,^ ,„ VnsA^nsT-^ o r dr 

[^T^r 3 Jo 



X I duj uj 
lo 



fl~U(r) 

"+5 5(-cc>) / d^G„+3(A*) 
Jo 



6.17 X IO^^^Lq nn{RTY T Jnpj.Tlf'^ 



'30 



(55) 



Values for the parameter 



1 



X 



n+3 



r^dr I dx- 
3 7o io 1 + 

X / d^Gn+sip) 
Jo 



(56) 



are shown Table iHIl 

The thermal cooling time scale of an isolated neu- 
tron star is of order 10^-10^ yr. The neutron star PSR 
J0953-)-0755 has an estimated age of 17 x 10^ yr, and 
HST observations indicate that it has a total luminos- 
ity of roughly 10~^Lq Comparing this to the total 
energy transferred to the star from KK graviton decay, 
Eq. H55|l . yields an upper limit on R. In Tables IVl and I VII 
we have tabulated limits on R and M4+„ derived from 
the maximum energy transfer per unit time to the neu- 
tron star Lniax = 1O~^L0. As we discussed in Ref. [T^ . 
this limit is far stronger than other astrophysical limits 
on R. For n < 5 it is stronger than any current labora- 
tory limits. 
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TABLE V: Upper limits on the compactification radius R (in meters) from our arguments. 

n 1234567 
Neutrino Signal 

SN 1987A 4.9 x 10^ 0.96 x 10"^ 1.14 x 10"^ 3.82 x 10~" 4.85 x 10"^^ 1.21 x 10"^^ 4.42 x 10"^^ 
EGRET 7-ray limits 

All cosmic SNe 4.9 0.96 x 10"^ 2.47 x 10'^° 1.21 x 10"" 1.93 x 10"^^ 5.6 x 10"" 2.29 x 10"" 

Gas A 4.1 X 10^ 3.61 x 10"^ 3.95 x 10""' 1.37 x 10"" 1.86 x 10"" 5.0 x 10"" 1.94 x 10"" 

Dv^T^occo^°,1^c;^. \ 1-935 2.63 X 10"** 7.17 x 10"" 3.91 x 10"" 6.99 x 10"" 2.24 x 10"" 9.98 x 10"" 
KX J 185635-3754 J 

Neutron-star excess heat 

PSR J0952+0755 4.44 x 10 



TABLE VL Lower limits on the fundamental energy scale M4+„ (in TeV) corresponding to the _R- limits of Table [V] Multiply 
limits on Mi+n by the factor M/Mn+4, given in the second row to obtain limits on the parameter M used in our previous 
papers [Hill. 



n 


1 


2 


3 


4 




5 


6 








7 


M/M„+4 


2.32 


2.98 


3.46 


3.82 




4.10 


4.32 








4.51 


Neutrino Signal 
























SN 1987A 


7.4 X 10^ 


8.9 


0.66 


1.18 


X 10"^ 


3.5 X 10"^ 


1.44 


X 


10" 


-2 


7.2 X 10"^ 


EGRET 7-ray limits 
























All cosmic SNe 


3.4 X 10^ 


28. 


1.65 


2.54 


X 10"^ 


6.8 X 10"2 


2.56 


X 


10^ 


-2 


1.21 X 10"^ 


Gas A 


7.7 X 10^ 


14.5 


1.24 


2.34 


X 10"^ 


7.0 X 10"^ 


2.80 


X 


10" 


-2 


1.37 X 10"2 


PSR J0953-H0755 \ 


4.58 X 10^ 


54.0 


3.46 


0.54 




0.14 


5.10 


X 


10" 


-2 


2.30 X 10"^ 


RX J185635-3754 / 








Neutron-star excess heat 
























PSR J0952-H0755 


1.61 X lO'^ 


7.01 X 10^ 


25.5 


2.77 




0.57 


0.17 








6.84 X 10"^ 



The predicted NS heating rate of Eq. (|55|) scales with 
the SN core temperature as T^'^/^, slightly steeper than 
for the previous arguments. The T-sensitivity of the lim- 
its on R and M^+n is almost the same as before. 



2. Reabsorption in the neutron star 

KK gravitons in the trapped cloud will be inside the 
NS on part of their trajectory. Assuming, as before, that 
all gravitons move on radial orbits we have calculated the 
fraction of time an average graviton spends inside the NS, 
^inside- While the graviton is inside the neutron star it 
can potentially be reabsorbed by the nuclear medium. 
Using Eq. pOfl we find the following absorption rates for 
average gravitons 



1 reabsorption ^71 P3 -^30 ' \^') 

where values of r„ are tabulated in Table |TV| Com- 
paring this reabsorption rate to the age of the two NSs 
RX J185634-3754 (1.6 x lO^^ s) and PSR J0953-t-0755 
(5.4 x lO^** s) we see that reabsorption in the medium is 
at most a moderate effect. 



V. CONCLUSIONS 

We have systematically revisited the constraints on KK 
graviton emission by SN cores and neutron stars that 
have been discussed in the literature. We have paid close 
attention to the scaling with the number n of extra di- 
mensions and to the conventions concerning the relation- 
ship between the radius R of the large extra dimensions 
and the fundamental energy scale, the A + n dimensional 
effective Planck mass M4+„ . Our limits are summarized 
in Tables ^ and IVII For n = 2 and 3 the limits based 
on EGRET observations are slightly less restrictive than 
stated in our previous papers, largely because we use a 
somewhat smaller emission rate from a SN core, and be- 
cause the estimated distance to RX J185635-3754 has 
increased by a factor of 2. However, the bound from 
neutron star excess heat is in fact slightly more restric- 
tive than stated in our previous paper. The reason is 
that the average distance of a graviton from the neutron 
star center is somewhat lower in the present calculation, 
compared to what was used in our previous work. 

The main theoretical new ingredient of our work is a 
new approach to calculating the emission rate. For low- 
energy gravitational bremsstrahlung the emission rate 
can be calculated directly from the classical formula in 
conjunction with the universal detailed-balancing prop- 
erty of the medium's dynamical structure function. This 
calculation agrees with previous works to the specified 
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order of precision, but is vastly simpler and clearly il- 
luminates the nature of the used approximation. The 
KK graviton emission rate by a SN core can not be cal- 
culated to better than a factor of a few. However, this 
uncertainty hardly affects the limits, especially for large 
n, because the limiting value of R is proportional to the 
n-th root of the limiting emission rate. 

Another benefit of our derivation is that it produces 
directly the absorption rate. It is found to be so small 
that the gravitationally trapped KK gravitons are not 
reabsorbed too quickly to invalidate our previous argu- 
ments. 
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Erratum "Supernova and neutron-star limits on large extra dimensions reexamined" 

(Dated: 13 October 2003) 



TABLE V: Upper limits on the compactification 


radius R (in meters) from 


our arguments. 






n 


12 3 


4 


5 




6 


7 




Neutrino Signal 
















SN 1987A 


4.9 X 10^ 0.96 X 10"'' 1.14 x 10"^ 


3.82 X 10~" 


4.85 X 10" 


12 


1.21 X 10"^^ 


4.42 X 10" 


13 


EGRET 7-ray limits 
















All cosmic SNe 


4.9 0.96 X 10"'^ 2.47 x 10""' 


1.21 X 10"" 


1.93 X 10" 


12 


5.6 X 10"^^ 


2.29 X 10" 


13 


Gas A 


4.1 x 10^ 3.61 X 10"'^ 3.95 x 10""' 


1.37 X 10"" 


1.86 X 10" 


12 


5.0 X 10"^^ 


1.94 X 10" 




PSR J0953+0755 \ 
RX J185635-3754 / 


7.36 5.13 X 10"^ 1.12 x 10^° 


5.46 X 10"^^ 


9.13 X 10 


-13 


2.80 X 10"^^ 


1.21 X 10 


-13 


Neutron-star excess heat 
















PSR J0952+0755 


4.44 x 10"^ 1.55 X 10"^° 2.58 x 10"^^ 


3.36 X 10"^^ 


9.95 X 10" 


14 


4.41 X 10"" 


2.46 X 10" 


14 


TABLE VL Lower limits on the fundamental energy scale M4+„ (in TeV) correspondini 


y to the 7?-Umits of Table El 




n 


12 3 


4 


5 




6 


7 






2.32 2.98 3.46 


3.82 


4.10 




4.32 


4.51 




Neutrino Signal 
















SN 1987A 


7.4 X 10^ 8.9 0.66 


1.18 X 10"^ 


3.5 X 10"^ 




1.44 X 10"^ 


7.2 X 10" 


3 


EGRET 7-ray limits 
















All cosmic SNe 


3.4 X 10=* 28. 1.65 


2.54 X 10"^ 


6.8 X 10"^ 




2.56 X 10"^ 


1.21 X 10 


-2 


Gas A 


7.7 X 10^ 14.5 1.24 


2.34 X 10"^ 


7.0 X 10"^ 




2.80 X 10"^ 


1.37 X 10 


-2 


PSR J0953-H0755 \ 
RX J185635-3754 J 


2.93 X 10^ 38.6 2.65 


0.43 


0.116 




4.31 X 10^ 


1.98 X 10"^ 


Neutron-star excess heat 
















PSR J0952+0755 


1.61 X 10^ 7.01 X 10^ 25.5 


2.77 


0.57 




0.17 


6.84 X 10 


-2 



In our paper "Supernova limits on large extra dimen- 
sions reexamined" [Phys. Rev. D 67, 125008 (2003)] a 
small error has occurred. We have mistakenly used the 
old distance of 60 pc rather than the new value of 117 pc 
for the nearby neutron star RX J185635-3754. There- 
fore, the corresponding limits on R are degraded by a 
factor (117/60)^/" while those on M4+„ need a factor 
(60/117)2/(2+«). Below we reproduce Tables V and VI 
with the corrected limits in boldface. We thank G. Sigl 



for drawing our attention to this error. 

Note that these limits can be significantly improved if 
one uses all neutron stars in the galactic bulge as a source 
rather than only the single nearby stars RX J185635- 
3754 and PSR J0953-h0755. The loss caused by the 
greater distance is significantly overcompensated by the 
large number of neutron stars. [M. Casse, J. Paul, 
G. Bertone and G. Sigl, "Gamma rays from the galac- 
tic bulge and large extra dimensions", hep-ph /0309 1 73( ] 



